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Highly ionized collimated outflow from HE 0238-1904 
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ABSTRACT 

We present a detailed analysis of a highly ionized, multiphased and collimated outflowing gas 
detected through Ov, Ovi, Neviii and Mgx absorption associated with the QSO HE 0238- 
1904 (Zem- 0.629). Based on the similarities in the absorption line profiles and estimated 
covering fractions, we find that the O vi and Ne vm absorption trace the same phase of the ab- 
sorbing gas. Simple photoionization models can reproduce the observed A^(Ne viii), A^(0 vi) 
and N{Mg x) from a single phase whereas the low ionization species (e.g. N iii, N iv, O iv) orig- 
inate from a different phase. The measured A^(Ne viii)/N{0 vi) ratio is found to be remarkably 
similar (within a factor of ~ 2) in several individual absorption components kinematically 
spread over ~ 1800 km s ' . Under photoionization this requires a fine tuning between hydro- 
gen density (ne) and the distance of the absorbing gas from the QSO. Alternatively this can 
also be explained by collisional ionization in hot gas with T > 10^^ K. Long-term stability 
favors the absorbing gas being located outside the broad line region (BLR). We speculate that 
the collimated flow of such a hot gas could possibly be triggered by the radio jet interaction. 
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1 INTRODUCTION 

Large scale outflows from AGNs play a vital role in regulating star 
formation in galaxies and the growth of the super-massi ve black 
holes at their centers dSilk & R ees 1998; King 2003; B ower et al.l 
I uOOq) . Hence detecting different forms of outflows is essential to 
■ understand the AGN feedback. In the QSO spectrum, outflows can 
manifest as associated narrow absorption lines (NALs) or broad ab- 
sorption lines (BALs) in the ultraviolet (UV) and as high ionization 
absorption lines and edges (i.e. warm absorbers, WAs) in the soft X- 
ray waveleng ths. BALQSO s comprise of ~ 40% of the total QSO 
population ( JDai et alJuOOSO while WA s are detected in the X -ray 
spectru m of ~ 50% of Seyfer t galaxies JCrenshaw et al.ll2003l) and 
QSOs jPiconcelli et alj|2005h . BALs & NALs are predominanfly 
detected through species with ionization potential (IP) < 100 eV 
(e.g. Civ, Siiv, Nv etc.), WAs, on the contrary, are identified by 
species with IP > 0.5 keV (e.g. O vn, O vm etc.). To understand the 
ionization structure of the outflowing gas entirely it is very impor- 
tant to detect species with intermediate ionization potentials such 
as Ovi (138 eV), Neviii (239 eV), Mgx (367 eV) etc. In partic- 
ular the resonant lines of O mAAlOSl, 1037 Ne vmAAUO, 780 and 
Mg xAA609, 624 are well suited for probing the intermediate phys- 
ical conditions of the outflowing gas. However the detections of 
such species in outflows have been very rare till date. 

There are three confirmed and two tentative detections of asso- 
ciated Ne vm (most of them show ing radio emiss i on) in the form of 
narrow absorption lines [UM 675. iHamann et al.l(ll995h ; 3C 288.1, 



iHamann et all ( l200d). J2233-6 6. [Petitiean & Sriananj (Il999l); 
HE 0226-41 lOr TGangulv et al] ( l2006h and 3C48. lOuptaetalJ 
( i2005') 1 and t hree QSOs i n the form of BAL absorption 
[Q 02 26-1024, iKorista et alj (Il992h: SBS 15 42+541, iTelferet al.l 
i ll 9981) and PG 0946+301. lAravetal.1 ( Il999l) 1. While multiphase 
photoionization models ar e generally used to explain these obser- 
vations (see for example, ' HamannI 119971) . the role of collisional 
ionization is not adequately explored. Here we report the detec- 
tion of a very strong associated Ne vm (and O vi) absorption in the 
high signal-to-noise (S/N) HST/COS and FUSE spectra of HE 
0238-1904 which has spectral energy distribution (SED) typical 
of a non-BALQSO. The system is at an ejection velocity of ~ 4500 
km s"' with a kinematic spread of ~2500 km s"'. Along with the 
Nevm and Ovi, we also detect absorption from Nm, Niv, Oiv, 
O v and Mg x. Throughout this paper we use cosmology with Q.^ 
= 0.27, Ha = 0.73 and /fo = 71 km s^'Mpc"'. 



2 OBSERVATIONS OF HE 0238-1904 

We use HS T/COS and FUSE FUV and NUV spectroscopic obser- 
vations of HE 0238-1904 that are publicly available in the MAST 
archivqj. The FUSE observations were carried out during Decem- 
ber 2000 (PI, Moos), July 2003 and October 2004 (PI, Howk). 



http ;//archive . stsci.edu/ 



S. Muzahid et al. 



We have downloaded a total of 86 frames calibrated using CAL- 
FUSE (v3.2.1) pipeline. The combined spectrum has a resolution 
of i? ~ 10, 000 and S/ N ratio of 3-10 per pixel over the wave length 
range 920-1180 A JMoos et al . 2000; Sahn ow et alj|2000l) . The 
archived COS observations are from the COS GTO Program 1 1541 
(PI, Green). These observations consist of G130M and G160M 
FUV grating integrations at medium resolution of R ~ 20, 000 and 
S/N ~ 20-30 per pixe l in the wave length range of 1134-1796 A 
JOsterman et a l.ll2011t i Green et al.l 12012). These spectra and the 
low-dispersion (R ~ 3000) NUV spectra covering the wavelength 
range of 1650-3200 A obtained with the G230L COS grating from 
the archive were extracted using the CalCOS v2. 15.4 pipeline. Con- 
tinuum normalization was done by fitting the regions free of ab- 
sorption lines with a smooth lower order polynomial. 

The optical spectrum of HE 0238-1904 was obtained with 
the 2-m telescope at lUCAA Girawali Observatory (IGOJl on 7-8 
December 2011. The long slit spectra (3x30 min exposures), cov- 
ering the wavelength range 3000 - 9000 A, were obtained with the 
GR5 grism of the lUCAA Faint Object Spectrograph (IFOSC) us- 
ing a slit width of 1.5 arcsec. The raw CCD frames were cleaned 
and the spectra were extracted using standard IRAF[j procedures. 
We use the method explained in Vivek et al. (2009) for removing 
the fringing at A > 7000 A. The wavelength and flux calibrations 
were performed using Helium-Neon lamps and standard star spec- 
trum respectively. The final co-added spectrum (in the heliocentric 
frame) has a resolution of R ~ 300 and a S/N ~ 20 - 40 per pixel. 



2.1 Spectral Energy Distribution and QSO parameters 

From simultaneous Gaussian fits to Hj8, O m AA 4960, 5008 emis- 
sion lines we get the emission redshift, Zem = 0.629 ± 0.002. 
The rest-frame UV/optical Spectral Energy Distribution (SED) of 
HE 0238-1904 is shown in Fig. \T\ The (green) star is the X- 
ray flux calcula ted from XMM-Newt on slew survey data assum- 
ing fy oc V •" l lSambrunaetal.lll999l) between 0.2 to 2 keV and 
Galactic A^(H i) = 3 x 10^°cm~^. Consi stency of photometric points 
(purple triangles) of Ojha et al. (2009) with our IGO spectrum (in 
red) implies lack of strong flux variability over 4.2 yrs in QSO rest 
frame. This is also supported by the similar UV flux measurements 
made by COS in December 2009 and STIS G140L in July 2002. 
In the absence of strong variability, multiple epoch data can be 
combined to get the QSO's SED. Fig. [T] shows that the source is 
intrinsically bringhter in UV tha n the HST compos ite spectrum of 
radio-loud QSOs constructed by ( iTelfer et al.l2002h . 

We measure X-ray-to-optical spectral index, a^x = -1.60 
which is consistent with a^y = -1.55 we get from the measure d 
optical luminosity L2500 using equation 8 of Stalin et alj ( I2OI0I) . 
This indicates that HE 0238-1904 is more like a typical non-BAL 
QSO. The fuU width at half maxima (FW HM) of the Hff lin e is 
~ 7200 km s"'. Using the prescription of iBentz et al.l ( 120090 we 
measure the size of the BLR to be ^^blr ~ 1.1 x 10'* cm. This 
together with th e FWHM of H/? gives the black hole mass Mbh 
= 2.4x10'" Mo dOnken et aLll20q41 and a Schwar zschild radius of 
i^sch ~ 7.0 X 10'^ cm. Following iHaU et alj J201lh . we find the di- 
ameter of the disc within which 90% of the 2700 A continuum is 
emitted D2700 ~ 3.2 x 10'' cm. 



^ http://www.iucaa.emet.in/~itp/ 

^ IRAF is distributed by National Optical Astronomy Observatories, which 
are operated by the Association of Universities for Research in Astronomy, 
Inc., under cooperative agreement with the National Science Foundation. 
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Figure I. Rest-frame SED of HE 0238-1904. The green star is from XMM 
slew survey data. The spectnim in cyan, blue, magenta and red are obtained 
from FUSE, COS(G130M, G16()M), COS(G230L), an d IGO respectively . 
The purple triangles are the photometric points from lOiha et al.l J20O9l) . 
The HST composite spectrum for radio-loud QSOs is plotted in light gray 
with arbitraiy normalization. Inset shows Gaussian fits to H/J, and Oiii 
/(/i4960,5008 lines which are used to estimate the emission redshift. 



3 ANALYSIS OF THE Ne vm ABSORBER 

In Fig.|2]we show absorption profiles of different species as a func- 
tion of outflow velocity with respect to the QSO (Zem = 0.629). The 
Neviii absoiption is detected in seven components spread over a 
~ 2500 km s"' velocity interval. The separate components are la- 
beled in the figure. Apart from the two weak Ne vm components 
(1 & 7), all other components show Ovi absorption with profiles 
very similar to that of Ne vm. The high ionization Mg xA624 and 
Ov/i629 lines are clearly detected in the FUSE data but without 
blending only in components 3 & 5. Due to poor spectral S/N the 
presence of stronger Mg x/1609 can be confirmed only for compo- 
nent 3. Contamination from Galactic H2 lines has made the de- 
tection of Mgx and Ov ambiguous in all other components. At 
the velocities of components 3 & 4, we also find weak absorption 
consistent with Nm/i989 and NiV/1765. The stronger N in/1685 line 
(fA ~ 1.7 times) is a non-detection in the FUSE spectrum, and 
therefore we report the A'(Nm) measurements for these compo- 
nents as upper limits. Very weak Ly/? and Lyy absorption are seen 
in the COS data for component 3. But we do not find any Lyo- ab- 
sorption in the low dispersion COS G230L spectrum in any of these 
Nevm components. 

The apparent column density profiles of Ne vm and O vi 
doublets, shown in Fig. [3] clearly indicates that the absorber is 
only partially covering the background emitting source. The es- 
timated covering fractions (/c ) using the method described in 
ISrianand & ShankaranaravananI (119991) are also shown in Fig. |3] 
The covering fractions estimated for the Ne vm in components 2-7 
is ~0.6. The Ovi absorption shows f^ ~ 0.5 in components 2 & 3 
and fc ~ 0.6 in components 4-6. The weak Nevm absorption in 
component 1 is consistent with a complete coverage of the back- 
ground source. It is interesting to note that Ne vm and O vi ions 
have very similar covering fractions for all the components. This 
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Figure 2. Absorption profiles of difl'erent species for the associated absorber shown with respect to the ejection velocity from the QSO. The Lyff is from low 
dispersion COS (G230L) spectrum whereas Mg x and O v are from FUSE. In total there are seven Ne vin components identified by the vertical dotted lines. 
The best fitting Voigt profiles after connecting for partial coverage are overplotted as smooth red curves. Absorption lines unrelated to this system are marked 
by the shaded regions. The blue smooth curve in the Mg x panel is the synthetic profiles with column densities as predicted from the photoionization model 
(see text). The vertical solid lines in the Mg x and O v panels show the positions of Galactic H2 lines. 



together with similar profile shapes suggest that the two ions are 
possibly tracing the same phase of the absorbing gas. The cover- 
ing fraction of f^ ~ 0.1 for Hi determined using Lyy6 & Lyy of 
component 3 is significantly smaller implying the kinematic coin- 
cide nce of multiple g as phases with different projected area (see 
also lXelfer et al.lll998l) . The results of Voigt profile fitting analysis 
after taking into account these effects of partial coverage is given in 
Table [T] In case of non-detections we calculate 3cr upper limits on 
the column densities using the rms error in the unabsorbed contin- 
uum at the expected position and by using the b value as measured 
for Ne vm (or O vi) with the appropriate f,.. 

The most intriguing fact about this system is that in the 
five components where we have measurements of A'(Ne vm) and 
A'(0 vi), we find a surprisingly constant value of 0.46±0. 16 for the 
log A'(Ne viii)/N{0 vi) ratio (see panel-D of Fig.[3}. To understand 
the physical conditions of the absorbing gas we ran several pho- 
toionization models with Cloudy v(07.02) (Ferland et al. 1998). In 
these models we assumed (a) the ionizing spectrum to be a com- 
bination of black body (with T ~ 1.5 x 10^ K) and power law 
with a-, = -0.7, a,,,, = -0.5 and q-o, = -1.6 as observed for 
HE 0238-1904 (see Fig.[T), (b) the gas is an optically thin plane 
parallel slab and (c) the re lative abundances of heavy elements are 
similar to the solar values JAsplund et al.ll2009n . 

In Fig. |4] various model predicted column density ratios are 
plotted against ionization parameter. In the case of components 
3 & 5, the observed A^(Nevm)/A'(Ovi), A^(Mg x)/A'(0 vi) and 
Af(Mgx)/A'(Neviii), within their 10% uncertainty, are well repro- 
duced for log U ~ 1.0 and 1.1 respectively. However, ratios involv- 
ing low ions (i.e. Oiv, Niv, O v) require lower values of log U in- 



dicating the presence of multiple gas phases. By using /^ = 0.6, the 
covering fraction determined for Ne vin, we obtain an upper limit 
of log A'(Hi) = 14.36 for the Hi column density associated with 
the high ionization gas. This gives a lower limit on metallicity of 
log (Z/Zq) > -0.8 for component 3. For all the components, the 
observed range of ^(Ne viii)/A'(0 vi) ratio is consistent with 0.95 
<logU< 1.30. 



Because of poor S/N of the FUSE data and the Galactic H? 
contamination we could only check the consistency of Mg x pro- 
files predicted by the models for other components. For each com- 
ponent we calculate the model predicted A'(Mgx) for log U re- 
quired to produce the observed A'(Ne vin)/A'(0 vi). The profiles are 
then generated assuming h and /^ as those measured for Ne vm. 
The predicted Mg x profile is consistent with the observed spectra 
within measurement uncertainties (see Fig. |2) suggesting similar 
value of log U for all these components spread over ~1800 km s"'. 



Nitrogen is reported to be over abundant compared to oxygen 
in associated absorbers fa voring a rapid enrichment scenario in the 
central regions of OSOs ([Hamann & Ferlandlll992l : iKorista et al.l 
1 19961 : |Petitiean&Srianandlll999l) . The best fit photoionization 
model predicts log Af(Nv) ~ 13.5. Even in a low dispersion spec- 
trum such a line should be detectable at ~ 4(T level. The fact that 
we do not detect N v line could be because (a) N is not over abun- 
dant compared to O and/or (b) the covering fraction of N v may be 
much less as it falls on top of Lyce-FN v emission lines. High reso- 
lution data is needed to confirm the above mentioned possibilities. 
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Figure 3. Panel (A) and (B) : appai'ent column density profiles of Ne vni 
and Ovi doublets respectively. The stronger and weaker transitions of the 
doublets are shown in red and blue respectively. Panel (C) : distiibutions of 
fc for O VI (in red) and Ne viii (in blue). Panel (D) : ^(Ne viii)/A'(0 vi) mea- 
sured in individual components labeled by the numbers. The mean and Icr 
scatter of the ratio are shown in horizontal dotted fines. The contamination 
are shown by the shaded regions. 



4 SUMMARY & CONCLUSIONS 

We report the detection of highly ionized outflowing gas through 
associated absorption from O vi, Ne viii and Mg x in the UV spec- 
tra of HE 0238-1904 which is presumably a non-BALQSO. The 
high S IN COS spectrum has allowed us to determine the cover- 
ing fraction {f^) and the multi-phase physical conditions in the ab- 
sorber. The similarity in the absorption profiles and fc suggest that 
O VI and Ne viii absorption possibly trace the same phase of the ab- 
sorbing gas. Near constancy of /^ for Ne viii and O vi in the different 
components spread over ~1800 km s"' indicates that the flow could 
be collimated. 

Under photoionization equilibrium conditions, the 
log A'(Ne viii)/A'(0 vi) = 0.46+0.16 ratio in the different compo- 
nents is consistent with 0.95 < log U < 1.30. This means that in 
spite of the absorbing gas being distributed over ~1800 km s"', the 
product nnr^ (where r^ is the distance of the gas from the QSO) 
is nearly a constant. Getting such constraint naturally for different 
components could be an issue for the photoionization model. 
Models of mass conserving shells expanding with a velocity Vj 
predict nnrjv, to be a constant. However, this is not the case in 
the present system as we do not find any clear trend between 
the ejection velocity and the ionization parameter. Alternatively, 
coUisional ionization in gas with temperatures in the range 
5.7 < log r(K) < 6.1 can recover the observed ^(Ne viii)/A'(0 vi) 
ratio as shown in the bottom panel of Fig |4] In this case we 
expect A'(Mg x)/A'(Ne viii) or A'(Mg x)/N{0 vi) to be very different 
between the different components, which cannot be ruled out 
with the existing FUSE data on Mgx. We note here that even 
in the case of photoionization the equilibrium gas temperature is 
~ 1.2xlO^KforlogU~ 1.0. 

From the observed flux at the Lyman limit (1.4x10"''* erg 
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Table 1. Partial coverage corrected Voigt profile fit parameters. 

Species Vrei(kms"') ' Up fo(kms"') log A'(cm"'^) J? 

Neviii -5767 ± 9 (T) 143.8 ± 12.7 14.22 ± 0.03 LO 

Hi .^ 143.8 <13.71 1.0 

Ovi -4743 ± 20 (2) 127.1 ±20.2 14.67 ± 0.09 OS 

Neviii 127.1 ±0.0 15.14 ±0.02 0.6 

Hi 127.1 < 13.70 1.0 

Ovi -4539 ± 6 (3) 98.4 ± 5.7 15.11 ±0.03 ois 

Neviii 96.7 ±1.9 15.45 ± 0.01 0.6 

Oiv 98.4 < 13.44 1.0 

Ov 98.4 14.45 0.5 

Nm 98.4 < 14.29 0.5 

Nrv 98.4 13.61 ± 0.04 0.5 

Mgx 98.4 15.32 0.6 

Hi ._^ 984 14.36(15.75) 0.6(0.1) 

Ne VIII -3678 ± 4 (4) 94.9 ± 3.5 15.10 ±0.02 (16 

Ovi 96.8 ±3.9 14.84 ± 0.02 0.6 

Oiv 96.8 14.39 0.6 

Nm 94.9 13.79 ±0.13 0.6 

Niv 94.9 13.55 ±0.04 0.6 

Hi .^ 94^9 < 14.11 1.0 

Neviii -3491 ±2 (5) 63.8 ± 2.9 15.11 ±0.02 (16 

Ovi 68.3 ±5.4 14.55 ±0.03 0.6 

Ov 68.3 14.03 0.6 

Mgx 68.3 15.11 0.6 

i 

Notes - ' Relative velocity with respect to Zcm- ^ Component ID as in 
Fig. |2] ' covering fraction used for the fit. Rest fra me wavelengths and 
oscillator strengths are taken from lVemer et al. I (119941 1. 



cm - s 'A ') in the QSO rest frame we get a relation log (unr^) = 
45.8 - log U. This gives log (uHr,-) = 44.8 for log U = 1.0. If we 
assume that the absorbing gas is well within the BLR (r;. < Rblr) 
we get Uh > 5 X 10* cm"^. The cloud thickness along the line of 
sight is, Ard < 2 x 10'' cm, if we use A?(Hi) ~ lO''* cm"^ and 
neutral fraction of /jjj ~ 10"^ as found for the photoionization 
model at log U = 1.0. The fact that the cloud is covering ~ 60% 
of the background emitting source gives the transverse size of the 
cloud as r^i ~ DxuMyjTd'^ = 1.2x10'^ cm~^. This is ~ 5 orders 
of magnitudes larger than the thickness along the line of sight (i.e. 
Arci/rci ~ lO""") resembling a sheet like geometry. 

If we assume spherical geometry for the absorbing gas then 
the observed f^ suggests that the radius of the cloud is ~ 1.2 x 10" 
cm. This gives Uh ~ iC cm"' and r^ ~ 90 pc for log U = 
1.0. Therefore the absorbing region, if nearly spherical, will be 
co-spatial with the narrow emission line region (NLR) gas. The 
sound crossing time is ~ 4 days and ~ 700 years respectively, 
for the two scenarios discussed above. Purely from stability point 
of view, the absorbing gas is likely to be outside the BLR (see 
[paucher-Giguere et al.l201 ll) .The NRAO VLA Sky Survey (NVSS) 
image shows possible radio emission associated with the QSO. We 
speculate that the outflow triggered by the radio jet could possibly 
explain the collimated hot outflowing gas outside the BLR as seen 
in the case of 3C48 (Gupta et al. 2005). 

The gas phase metallicity is found to be, Z > 0.2 Z© with- 
out a clear evidence for an enhanced N abundance. However, high 
resolution spectra covering C iv and N v is needed to perform de- 
tailed abundance analysis. High densities (uh) inferred for the sys- 
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Figure 4. Top: Results of photoionization calculation. In each sub-panel the 
solid curve shows the Cloudy predicted column density ratios against log 
U. In the bottom most panel the measurements of WCNe viii)/N(0 vi) ratios 
in five components are shown by the stars together with their component 
ID. For other panels ratios are shown only for component 3 & 5. Bottom: 
Column density ratios versus tem perature under collisional ionization equi- 
libiium fOnat & Stemberg|2007b . The horizontal shaded region shows the 
range in the measured N(Ne vui)/N(0 vi) whereas the vertical shaded region 
shows con'esponding allowed range in temperature. Clearly within the per- 
missible range of temperature, ^(Mg x)/N(0 vi) and/or A'(Mg x)/A'{Ne vm) 
ratios can vary more than an order of magnitude. 



tern suggests that HE 0238-1904 is an ideal target for absorption 
line variability studies. However, we do not find any signature of 
strong variation in Ne viii and O vi equivalent width from 2002 on- 
wards from the HS T/STIS and COS spectra. As this source is X- 
ray bright detection of warm absorbers will allow us to get further 
insight into this hot outflowing gas. 
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